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SECTION  1 
INTRODUCTION 


This  report  is  the  first  in  a series  of  reports  to  appear  on 
analytic  studies  of  general  SGEMP  phenomena.  They  are  intended  to  investi- 
gate the  physical  phenomena  involved  in  the  generation  of  skin  currents, 
surface  charge  densities,  and  electromagnetic  fields,  and  in  coupling 
mechanisms  to  internal  circuitry.  It  is  hoped  to  provide  at  least  a 
qualitative  and  preferably  a quantitative  understanding  on  theoretical 
grounds,  to  provide  a physical  understanding  of  the  various  computer  code 
results  obtained  to  date,  and  to  help  in  predicting  expected  order  of  magni- 
tude effects  on  real  three-dimensional  targets. 

The  spirit  of  our  approach  is  to  adopt  Fermi's  philosophy  that 
"one  should  never  do  a (detailed)  calculation  until  one  already  knows  what 
the  answer  is."  We  wish  to  outline  the  basic  physics  of  SGEMP  phenomena  so 
that  we  have  a sound  theoretical  framework  in  terms  of  which  we  can  under- 
stand the  more  detailed  computer  results. 

The  present  report  deals  with  the  steady  state  electron  emission 
boundary  layer.  In  the  future  we  hope  to  study  the  time-dependent  transient 
build-up  of  the  layer. 

T*  we  think  of  the  chronological  sequence  of  events  occurring  in 
a typical  SUkMP  problem,  we  have  first  the  arriving  X rays  illuminating  the 
exposed  satellite  surfaces  and  ejecting  photoeloctrons.  These  electrons 
constitute  an  electric  current  which  is  the  basic  driver  of  Maxwell's 
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equations.  The  exposed  surface  is  left  positive  and  this  positive  charge 
flows  away  on  the  satellite  surface  producing  the  skin  currents. 

The  photoelectrons  above  th •:  positive  surface  produce  a strong 
noraal  electric  field.  At  relatively  high  fjluences,  large  numbers  of  low- 
energy  electrons  are  held  back  by  this  field  and  congregate  near  the  surface, 
producing  a boundary  layer  of  relatively  dense  charged  plasaa.  Of  the 
subsequently  eaitted  electrons,  only  the  more  energetic  ones  can  penetrate 
the  potential  of  this  layer  and  move  through  space  around  the  satellite. 

Hence  the  boundary  layer  greatly  affects  electrons  that  penetrate 
it  by  reducing  their  energy.  Also,  the  layer  constitutes  a dipole  layer 
producing  its  own  external  field.  In  addition,  the  structure  of  the  layer 
near  the  edges  of  exposed  surfaces  will  determine  how  the  electrons  flow 
around  the  edges,  and  may  help  in  understanding  the  surge  of  replacement 
current  around  the  sides. 

The  boundary  layer  structure  is  a function  of  the  X-ray  fluence, 
time  history,  and  energy  spectrum,  and  the  photoelectric  yield  of  the  target. 
Its  thickness  can  be  less  than  a millimeter  for  high  fluences  and  soft 
spectra,  and  larger  than  the  satellite  dimensions  for  hard  spectra  and  low 
fluences  so  that  it  hardly  exists.  This  wide  range  of  thicknesses  makes  it 
difficult  to  handle  by  a single  technique  in  SGEMP  codes  designed  for  a 
wide  range  of  fluences  and  spectra. 

In  this  report  we  study  the  structure  of  the  boundary  layer  in 
the  steady-state  approximation  and  in  one  dimension.  In  this  case  the 
problem  can  be  reduced  to  quadratures,  and  so  solved  completely.  Previous 
workers  have  investigated  various  aspects  of  this  problem1 »2 * 1 . Here  we 
consolidate  these  previous  results  from  a unified  point  of  view  and  present 
some  new  results  of  interest. 
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SECTION  2 

CHARACTERISTIC  DIMENSIONS 


Figure  1 shows  a sketch  of  the  relevant  geometry.  The  z direc 
tion  is  normal  to  the  surface.  The  X rays  are  incident  from 
above  in  the  normal  direction.  Electrons  are  ejected  upward  with  some 
angle  and  energy  distribution.  In  the  steady  state,  the  electron  number 
density  n(z),  electric  field  E(z),  and  potential  4>(z)  will  all  be  single 
valued  functions  of  z. 


Let  Wj  be  a characteristic  energy  of  the  emitted  electrons,  and 

(1) 


v j » V - wi/m  * 


he  a characteristic  velocity.  The  electron  number  density  at  the  surface, 
n(0),  which  includes  both  emitted  and  returning  electrons,  defines  a plasma 


frequency  u>  P ^ 


EZ 

Wp  y in 


n(0)  . 


(2) 


The  energy  w^  defines  the  characteristic  dimension  the  Debye 


length,  by 


i 


V? 


(3) 


u 


i(OJ  --p 

This  is  the  distance  over  which  n,  E,  and  vary  appreciably. 
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Figure  1.  Sketch  of  geometry. 

The  other  dimensional  unit  of  distance  is  e2/w. . But  this  is  not 

1 2 

a plasma  quantity  and  is  at  any  rate  exceedingly  small.  Since  e * 1.44  x 
10'10  keV  cm,  and  is  a few  keV,  e2/Wj  < 10”^  cm,  and  is  not  a dimension 
of  concern.  It  is  the  distance  of  closest  approach  of  two  electrons  with 
relative  energy  Wj. 


The  average  distance  between  electrons  is 


(4) 


If  electrons  are  ejected  at  a rate  rQ  (electrons/cm2/sec) , which 

is  the  product  of  the  material  yield  Y(elecs/cal)  and  the  X-ray  flux  dd»/ dt 
2 

(calories/cm  /sec), 


electrons 

2 

cm  sec 


(5) 


then  the  number  density  at  the  surface  will  be  on  the  order  of 
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n(0) 


(6) 


— — eiecs/cm^  , 

V1 

and  the  total  number  of  electrons  will  be  on  the  order  of 

/mV  0 2 

Ng  ~ n(0)^D~  j-  elec  s/cm  . (7) 

V 4ne 

The  surface  charge  density  a will  be 

a * eNg  y mv jTQ/4TT  esu/cm^  , (8) 

and  is,  curiously,  independent  of  the  electron  charge.  If  e were  not  a 
constant,  for  example  if  it  were  larger,  each  ejected  electron  would  leave 
behind  more  charge,  but  the  electrons  would  return  more  rapidly  to  the  surface, 
so  that  in  the  steady  state  a does  not  depend  on  e. 

Likewise  the  surface  electric  field,  ECO),  will  be  of  order 

ECO)  = 4tt o - V7rrmvp:“'  esu  , (9) 

1/2 

and  is  also  independent  of  e.  Both  a and  E(0)  increase  as  r^  and  there- 
for. with  the  square  root  of  the  fluence. 

The  potential  of  the  surface  relative  to  a distant  point  will 
be  of  order 

w 

*1  ~ E(0)JLd  - 4TTen (0)£p  = — . (10) 

Thus  the  characteristic  potential  depends  only  on  the  electron  energy 
spectrum,  and  is  independent  of  fluence. 

Tuc  dipole  moment  per  unit  area  of  the  boundary  layer  will  be  of 

order 


p 


oN  = a2,n~  — ~ 
s 0 D 4Tte 


csu-cm/ cm“ 


(11) 
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also  independent  of  fluence.  Thus  one  might  expect  that  the  quasi-static 
field  produced  by  the  boundary  layer's  dipole  moment  will  become  relatively 
less  and  less  important  as  fluence  increases  since  this  field  is  independent 
of  fluence  whereas  other  sources  of  dynamical  fields  increase  with  fluence. 

The  kinetic  energy  stored  in  the  moving  electrons  will  be  of  order 
1 2 / mVo  2 

K ~ w.N  ~ -T  mvi  / T ergs/cm  , (12) 

15  2 1 V 4i.e2 

and  the  energy  stored  in  the  electric  field  will  be  of  order 

2 

u ~ ML-  ~ K ergs/cm2  . (13) 

Later  we  show  that  general  principles  actually  imply  that  IJ  = 2K  exactly 
for  the  normal  component  of  motion.  Subsequent  detailed  calculations  in  this 
report  confirm  all  of  the  above  dimensional  arguments. 

TIME  CONSIDERATIONS 


The  characteristic  turn  around  time  for  an  electron  trajectory  will 
be  of  order 


1 


1 


a) 


I mv] 

2 

8tic  r 


(14) 


0 


.13 


We  estimate  tj  when  Wj  = 2 keV,  Y = 1.25  x 10  elec/cal  (2  keV  blackbody 
on  Aluminum),  and  the  X-ray  pulse  FWUM  is  3 shakes.  Then 

I iT1” 

1.4  / -nr-  nanosec  , (15) 


’I 


'(Hcal/cm  ) 


where  4>  is  the  X-ray  fluence,  and  we  have  approximated  d<|>/dt  in  liquation  5 

- 3 "> 

by  its  average  value  <p/2  FWHM.  Hence  for  fluenccs  above  10  ' cal/enT, 

tj  < 1 ns.  This  indicates  that  the  steady-state  solutions  discussed  in  this 


8 


I 


report  will  be  applicable  when  the  X-ray  flux  is  not  changing  appreciably 

-S'*  . i ■* 

over  times  of  order  1 ns  at  10  cal /cm"  or  of  order  0.1  ns  at  10  cal /cm". 


The  steady-state  solution  should  begin  to  become  a good  approxima- 
tion when  the  total  number  of  electrons  ejected 

t 

N(t ) = V f dt  elecs/Cu"  . 


(16) 


0 


is  larger  than  the  steady-state  number 


,‘(t>  - v & (t)  • 


(17) 


from  Equation  7.  Assuming  a triangular  time  pulse  with  a rise  time  tf. 
and  a full  pulse  width  T = 2 FW'IJM, 


d£ 

dt 


£ _ 

T t 


_t 

rise 


0 < t 5 t 


rise  ’ 


(18) 


for  the  rising  portion  of  the  pulse.  Using  18  and  equating  lb  and  1“  one 
obta i ns 


mv 


f = t 


rise 


2-e 


L _L 

: v: 


{ i9i 


or 


_ 1 x 


. :s 


V w . ( ke\  ) t . i shakes)  l l shakes) 
1 rise 


.15 


_ _ 1 / 3 


Using  Y = 1 C5  \ It)  ' elec  cal,  T = shakes,  t 

r l se 


t = 1 


■>( 

\ : i ca 1 /cm  ) / 


nanosec  , 


ns  . CO) 

1 shake,  w * 2 keY , 

Cl  I 


which  varies  from  10  to  1 ns  as  : increases  from  10  1 to  10  * c.»l/cm“. 
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SECTION  3 


THE  RELATION  BETWEEN  EMISSION  ANGULAR  DISTRIBUTION 
AND  NORMAL  ENERGY  DISTRIBUTION 


In  our  one-dimensional  problem  only  the  normal  component  of  motion 
is  affected.  Electrons  are  ejected  with  some  energy  and  angle  distribution 
d“n/dwd£l  (electrons/cm  sec  kcV  ster) , but  the  dynamics  influences  only  the 
z component  of  velocity  v„.  The  emission  energy  distribution  is 


dn 

dw 


_dji_ 

dwdft 


da 


electrons 

i ► 

cm“  sec  keV 


(22) 


with  the  integral  taken  over  ^ne  emission  hemisphere,  and  the  three-dimensional 
velocity  distribution  is 


m d"n 
v dwda 


electrons 
c»“  sec  (cm/sec)^ 


(23) 


The  z component  of  velocity  is  distributed  according  to 


where 


— r~ 

\ " v v' 


i:4) 


U'Sl 
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and  $ is  the  azimuth  angle.  We  now  assume  the  emission  spectrum  factorizes 
into  energy  and  angle  parts. 


si®s  ■ sr  ’ 


(26) 


where  9 is  the  polar  angle  from  the  normal  as  in  Figure  1.  Then  using 
dw  * mvjdvt,  valid  when  is  held  fixed, 


— a.  ii 

dn  _ C dw  dn  C . x. 


(27) 


where 


1 2 
w * -=■  m v 
Z 2 2 


(28) 


is  the  normal  component  of  energy.  We  now  assume  the  emission  is  a cosd 
distribution. 


g(3,<?)  = £ cosO  = ^ , 

which  is  close  to  the  experimental  facts  and  to  the  distribution  predicted 
by  the  code  QUICKE2.  In  this  case.  Equation  27  yields  for  the  distribution 
of  w , 


L = / electrons  \ 

/ w dw  I 2 . ) 

z J ‘cm  sec  keV ' 


dn  _ 
dw 


(29) 


The  following  several  examples  of  energy  distributions  dn/dw  with  a cost? 
angle  distribution  and  the  resulting  dn/dw,  are  instructive. 

1 . Monoenerget ic 


dn  , . 

as  = V f* 


(30) 
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(31) 


-^2-  « e(w,  - w ) , 
dw  w,  v 1 zJ  * 
z 1 

where  0(x)  is  the  unit  step  function  0 (x)  » 1 if  x > 0,  6(x)  ■ 0 if  x < 0. 
A monoenerget ic  cosG  distribution  produces  a constant  distribution. 

2.  Constant 


dn 

dw 


w)  , 


0 < w < w.  , 
z 1 


3.  Linear  Increase 


dn 

dw 


0 < w < Wj  , 


dn 

dw__ 


M 


(W1 


*0 


0 < < Wj 


(32) 


(33) 


(34) 


(35) 


A linearly  increasing  cos©  distribution  produces  a linearly  decreasing 

4m 

distribution. 


4.  Linear  Decrease 


dn 

dw 


w) 


0 < w < Wj  , 


dn 

dw 


w 


0 < w_  «-  Wj 


5.  Exponential 


dn 

dw 


-w/  w. 


(36) 


(37> 


(38) 
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(39) 


CO 


■ sfWV  • 


where  Ej(x)  is  the  exponential  integral.  This  case  is  of  interest  because 
actual  electron  energy  spectra  when  blackbody  photon  sources  are  incident 
are  very  nearly  exponential4.  Note  dn/dwz  diverges  logarithmically  as  w2  ■*  0. 

6.  Linear  Times  Exponential 


dn 

dw 


-W/Wj 

w e 


(40) 


dn 

dw 

z 


ro  -V“i 

=7* 


(41) 


Here  the  z "component"  of  energy  is  exponentially  distributed.  The  energy 
spectrum  40  is  essentially  an  exponential  with  the  low-energy  electrons 
deleted.  Hence,  comparing  results  for  cases  S and  6 above  will  indicate  the 
effect  of  the  low-energy  electrons.  For  example  we  shall  find  that  the 
slope  of  the  charge  density,  dp/dz,  at  z * 0 diverges  for  case  S but  is 
finite  for  case  6.  Equation  40  is  also  the  emission  energy  spectrum  for 
particles  of  a Maxwellian  gas  escaping  through  a small  hole. 


These  six  energy  spectra  are  sketched  in  Figure  2.  The  cosO 

distribution  always  enhances  the  low  energy  end  of  dn/dw„.  In  particular 

*• 

dn/dw^  diverges  logarithmically  as  w . -*  0 if  dn/dw  is  non-zero  as  w ■*  0. 
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Figure  2.  Normal  energy  distributions  from  several 
cose  distributions. 
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SECTION  4 

SOLUTION  FOR  BOLTZMANN  DISTRIBUTION  FUNCTION 


The  problea  aay  be  set  up  by  solving  for  the  scalar  potential  ♦ 
fro*  Poisson's  equation 


« - 4irp(z)  . 
3 z 


(42) 


and  for  the  particle  Motion  fro*  the  tine  independent  Boltz*ann  equation 
for  the  distribution  function  f 


where 


and 


v • fa  - J « • j£  - J B(»  g-  . 0 


E(z)  - - 5J-  . 


(45) 


PU)  " - 


d3vf  . 


(44) 


(45) 


By  dividing  43  by  -evzE(z),  and  using  E(z)dz  ■ -d$,  and  dw2  ■ *vzdv2>  we 
obtain 


3f  3f 

e5*  * to~z  " 0 ’ 


(46) 


the  general  solution  to  which  is  any  function  of  w,  - e4>, 

f = f(w,  - c<2>)  . (47) 
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The  correct  function  to  use  in  17  is  obtained  by  matching  the  boundary 
condition  on  f at  the  emission  surface,  obtained  as  follows. 


The  differential  rate  of  emitting  electrons  is  f(z  » 0)vz  and  is 
given  by  Equation  23, 


f(z  ■ 0)v 


3 2 

d n b m d n elocs 

z " .3  v dwdfi  2 , , ,3  * 

d v cm  sec  (cm/sec) 


(41) 


Again  assuming  a cos6  distribution,  we  have 


f(z  • 01  ■ jr  ac . cos6 

Z 


rrv 


m dn 
2 dw 


m 1 dn 
2ir  w dw  * 


(49) 


The  distribution  of  v is  obtained  by  integrating  out  the  v and  v dependence 

z x y 

to  the  one-dimensional  distribution  function 


so  that 


^(^’V  * j dvxdvy  f(z*^  * 


flO 


O.v,)  - j i 


dv  dv  f(C,v) 
x y ' • • 


(50) 


/dw  dn 
w dw 


(51) 


Since  Equation  16  also  holds  for  fj,  47  implies,  with  the  help  of  29, 


fiu-v 


, f dw  dn 
J w dw 


w - e<!> 
z 


dn.  lwI  - o») 
dw 


(52) 
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We  have  Multiplied  by  a factor  of  2 in  Equation  52  to  account  for  the  return- 
ing electrons  in  the  steady  state,  and  have  taken  $(z  • 0)  ■ 0. 


The  charge  density  45  is  then 


p(s)  * - e fj 

0 


* - e 


dw2  _dn_ 


(53) 

4 


• P(0)  • 

An  alternate  expression  can  be  obtained  with  jjJ-  by  using  the  first  line  of 
liquation  52  and  interchanging  the  order  of  integrations  on  v„  and  w, 

4* 

an 

PU)  * - 2c VI*  j Vw  ♦ e*  . (54) 

-eif 


These  arguisents  and  some  in  the  following  section  are  similar  to  those  of 
Higgins’ . 
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SECTION  5 

SOLUTION  OF  POISSON  EQUATION 


With  p given  by  S3  or  54  as  a function  of  <t>,  Poisson's  equation  is 


a2* 

— =-  ■ - 4trp{$) 

3Z*4 

i a / 1>4> 

" 2 W3z) 


(55) 


so  that 


E « - 3<J>/3z 

4>(z) 

■ ! - 81T 


I 


J p($')d$' 


♦ E(a) 


,|l/2 

I ’ 


(56) 


*(«) 


is  the  electric  field  as  ■»  function  of  G'(z).  Here  E(a)  is  an  integration 

constant,  the  field  at  any  point  z • a.  If  the  emission  energy  spectrum  has 

a maximum  energy  w„(max),  the  boundary  layer  will  have  a maximum  extent 

and  the  maximum  potential  will  be  $ » -w  (max)/c  = $U  )»  ’»nd  E ( : ) = 0. 

aid  x z max  max 

Thon  the  electric  field  *\11  be 


E(z) 


♦ 

max 


t 


(57) 


If  the  emission  energy  spectrum  n..  '^cvrons  of  arbitral  i 'v  high  energy,  then 

z = 00 , li(°°)  = 0,  and  liquation  57  will  still  hold  with  £'  - - Integrat- 

max  max 

ing  once  again, 
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2 ■ - 


♦ 

^ 1 

/ d*'  | 

- Sir  1 p($")d$" 

J 

•/  i 

0 < 

* 

-1/2 


max 


(58) 


The  Inverse  of  S8  gives  $(z)  and  completely  solves  the  one-dimensional 
steady-state  problem.  This  solution  was  given  by  Hale  for  several  spectra2. 


Equations  r.7  and  S3  give  for  the  electric  field  at  the  surface 

Q w2(max) 

E(0)2  - ftxVS/df  f ■ tS9) 

A J Z 2 

max  -e$* 

Interchanging  the  orders  of  integration  permits  the  4>*  integral  to  be  done. 
Let  Wj  be  a characteristic  energy  in  the  emission  spectrum.  Then  we  find 

wz(max) 

E(0)2  = 167rV35^  r0  j . (60) 

0 

The  integral  is  a dimensionless  function  of  the  emission  spectrum,  and 
hence  Equation  60  confirms  our  dimensional  arguments  in  Equations  9,  8,  and 
7,  for  E(0),  o * E(0)/4tt,  and  •■a/e. 

The  dipole  moment  per  unit  area  P(z)  contributed  by  all  electrons 
out  to  a distance  z is 

z 

P(z)  ■ J zp(z)dz  . (61) 

0 

Using  3E/3z  = 4irp  and  integrating  by  parts,  this  can  be  expressed  as 

P(3)  * ^ (*(Z)  ♦ zE(z)]  . (62) 

If  the  emission  spectrum  has  a maximum  energy  wz(max),  the  total  dipole 
moment  will  be 
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P - PU( 


(63) 


$ 

, b max 
maxJ  “ 4tt 


wz(max) 

4ie 


independent  of  fluence  as  in  Equation  11.  If  there  is  no  maximum  energy, 
P(z)  may  diverge  as  z ®.  This  occurs  for  cases  S and  6,  and  P(z)  is 
only  weakly  dependent  on  fluence  when  z >> 
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SECTION  6 


ENERGY  CONSIDERATIONS 


The  equations  of  motion  for  the  electrons  can  be  written  in  the 
form  of  a conservation  law  by  introducing  the  momentum  densities  tt  and  (j, 
and  momentum  flow  tensors  and  of  the  particles  and  fields*.  Then 
these  equations  read 


where 


Jr  (it.  + G.)  ♦ — (P.  . + T.  .)  * 0 , 
3t  v 1 i'  v ij  ly  ' 


fd3v 


\ ■ 

P . = / mv . v . 

ij  J 1 J 

G.  = (£  x t). 

1 47TC  1 


fd3v 


(64) 


1 r 1 


Tij  = 4?  ff  <E  * B >«ij  - <EiEj  + BiVJ  • 

In  steady  state  one  dimension,  the  time  derivative  and  x and  y derivatives 
in  64  vanish,  leaving 


Vz  (Pi3  * Ti3>  = 0 ' 


(65) 


so  that  P.,  + T._  is  independent  of  z.  In  the  present  problem,  E.  = E5 . _ , 

1 J 1 j 1 1 J 


B = 0,  so  the  i = 3 component  of  65  reads 

-> 

2 3 £•■ 

mv.  fd  v > ^ = constant 


/■ 


(66) 
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At  large  z,  both  terms  on  the  left  hand  side  of  66  separately  vanish,  so 
the  constant  is  zero.  The  energy  density  stored  in  the  electrostatic  field 
is 


and  the  density  of  z component  of  kinetic  energy,  k , is 

z 

1 fl  2 £.l 
k2  = J y n>vzfd  v , 


(67) 


(68) 


so  that  66  implies 

UE  " 2 \ ■ 


(69) 


The  electrostatic  energy  density  is  everywhere  twice  the  z kinetic  energy 
density.  With  a cos9  emission  distribution,  f is  independent  of  angles  so 
the  x and  y contributions  to  the  kinetic  energy,  k , and  k , are  each  equal 

x y 

to  k . The  total  kinetic  energy  density  is 


k 


■ f?m 


y 

2 


+ k 

z 

,.3 

f d v 


(70) 


and  Equation  69  implies 


(71) 


Integrating  on  z,  the  total  stored  electrostatic  energy 

'’max 

UE  = f uE(z)  dz  , 

0 


and  the  total  stored  kinetic  energy 
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are  related  by 


K 


max 

J k(z)  dz  , 
0 


1 


(72) 
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SECTION  7 


SEVERAL  EXAMPLES 


In  this  section  we  collect  some  formulas  and  graphs  for  four 
interesting  examples. 


Example  1 - Monoenergetic  Normal  Emission 

2 

Electrons  are  emitted  normally  at  r^  elecs/cm  /sec  all  with  energy 

I ” 

W1  " yEVl‘  This  carc  was  ^irst  discussed  by  Karzas  and  Latter1. 


The  number  density  at  : = 0 is 
n (0)  - 2r(J/v1  , 

and  the  Debve  length.  Equation  3,  is 


rvT 


The  maximum  distance  the  electrons  travel  is 


“max  3 

The  electric  field  at  the  surface  is 
E (0)  = 4 V Tinryv  j 


3 e: 

max 


(73) 


(74) 


(75) 


(7b) 
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The  surface  charge  density  is 


a - vWqV^tt  , 

and  the  total  number  of  electrons  is 


(77) 


_ / mr  v 

N = — a -A 


, p-  * 3 n(0)  2 /— e?-\ 

s e V Tte2  max\  cm2  / 


The  number  density  is 


(78) 


n(z)  = n(0)  (1  - 2/zmax) 


-2/3 


(79) 


and  diverges  as  : ■*  z . The  e^ectr^c  field  an<i  potential  are 


E(:)  * E(01  (1  - z/z  ) 
v ' v . v max' 


1/3 


♦«*)  •r'11  - ;/WJ/5  - *'  • 


(80) 

(81 


and  the  dipole  moment  is 

w 


P = -pi-  = 4 eN  = esu/cm  . 

4*e  4 s max 


(82) 


The  stored  kinetic  energy  is 


K * j w1Ns  . 


(83) 


and  the  stored  electrostatic  energy  is 


U,;  = 2 K , 


(84) 


since  there  is  no  x or  y motion. 


The  time  for  an  electron  to  reach  ; is 

max 


5 ; / mv, 

t > — Hi  = /-^~ 

max  'l  V 4~e“r, 


(85) 
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An  electron's  trajectory  is 


V(t)  - »,(1  - 

0 < t < t , 

max 

(86) 

z(t)  « z „ [1  - (1  - t/t  )3]  0 < t < t 

v max1  v max  J max 

(87) 

and 

v(z)  - Vj(l  - z/zmax)2/3 

(88) 

An  electron 

leaving  the  surface  sees  an 

electric  field  decreasing 

linearly 

with  time. 


The  number  density  n(z)/n(0)  and  electric  field  E(z)/E(0)  normalized 
to  their  value  at  z = 0 are  shown  as  curve  1 in  Figures  3 and  4.  The  normalized 
potential,  -e<t>/Wj»  is  shown  as  curve  1 in  Figure  5. 

Example  2 - cos9  angular  distribution,  monoenergctic  energy  distribution 
(case  1 of  Section  3).  This  is  the  same  as  a normal  emission, 
constant  spectrum,  and  was  discussed  by  Hale2.  Here 


dn  j-,  . 

a?*  V(w  - wi>  • 

and,  with  Vj  = V 2Wj/m,  we  find 


n(0)  = 


4r. 


1 


, . /-Vl_ 

D V 16ne“r0 

W ;VTin  • 


E (0)  = 4 \j  Jirmr^v  j / 3 


- 4 


max 


(89) 

(90) 

(91) 


(92) 
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Figure  3.  Normalized  electron  number  density  vs.  distance 
from  surface. 
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* 

s ipyr r Kip*  s to " n m-  117* S’-1 * re • 


.^■.v,L.k-->..  tteaxto  BiArt  jtfc^a  Jia^t  ^ rlfaWt: 


o • V^argVj/SJt  , 


»s  *1*  ji>(°)w . 


n(z)  - n(0)  (1  - z/z„„r  . 


E(z)  - ECO)  (1  - Z/2  J3  . 


max' 


♦“>  ■ T 1(1  ' - »)  • 


max' 


(93) 

(94) 

(95) 

(96) 

(97) 


P * 4^  * 7 eNszmax 
Kz  “ 7 wlNs  » 


(98) 

(99) 


UE  ‘ 2 Kz  “ f K » 


v(z)  - Vj(l  - z/z^) 


(100) 

(101) 


The  maximum  distance  an  electron  with  z component  of  energy  w , 

0 < wz  < w , will  go  is 


Wwz}  * 2aax 


(1  - (1  - wywj)174]  . 


(102) 


Trajectories  and  turn  around  times  involve  elliptic  functions,  but  in  the 
case  of  a particle  emitted  with  energy  - Wj/2,  we  find 


W"i  * 7 V * °-49  -ff 


(103) 


The  number  density,  electric  field,  and  potmtial  are  plotted  as 
curve  2 in  Figures  3,  4,  and  5. 
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i 

i 


Example  3 - cosO  angular  distribution,  linear  times  exponential  energy 

distribution  (case  6 of  Section  3).  In  this  example  the  inte- 
grals in  Equations  57  and  58  can  be  done  analytically. 


dn 

dw 


-w/w, 


w e 


We  find 


^2wj/m  . 


Define 


n(0)  * 2VT  ~ , 
V1 

rv: 

S - n £d  * 


(104) 

(105) 


which  is  the  Debye  length  at  the  surface  of  the  electrons  which  are  leaving 
the  surface.  Then 


E(0) 

* V 8Tr5^mrQvl 

2"i 

■*r- 

(106) 

0 * “ 

VmrQVj/Stnr  , 

(107) 

N = 
s 

“ * n(0U.  , 
e l 

(108) 

n(z) 

- n (0)  (1  ♦ z/tj)'2  , 

(109) 

E(z) 

= n (oi  (i  ♦ z/f.j)'1  , 

(HO) 
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(Ill) 


2wl 

*(z)  - — Jtn  ci  ♦ i/ly)  , 


Kz  “ f wlNs  • 


UE  “ 2Kz  "I  K * 


(112) 

(113) 


The  dipole  moment  per  unit  area  due  to  electrons  out  to  z is 

* as  [*"»  * «/*i>  - rnq-]  • 


(114) 


For  z » Jlj  this  depends  only  logarithmically  on  the  fluence  through  tj. 


1 2 


An  electron  emitted  with  velocity  v and  energy  w » mv  has 


velocity 


f*  w 

V(z)  - v[l  - 2 in(l  ♦ z/1 lL)]1/2  , 


(US) 


and  reaches  maximum  distance 


/ w/2wl  \ 

- - V ' 


z (v) 
maxv 


(116) 


in  a time 


tr=  w/2w 

t (v)  • erf(v/VTv.) 

maxv  y 2 v^  v 1' 


(117) 


where  erf(x)  is  the  error  function.  For  v ■ v^ , 

WV  * l-405  ^ • 


(118) 


Figures  3,  4,  and  5 show  n(z),  E(z),  and  4>(z)  in  curve  3. 
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Example  4 - cosd  angular  distribution,  exponential  energy  distribution 


(case  S of  Section  3) . This  distribution  corresponds  closely 


with  the  experimental  distributions  for  blackbody  photon 
The  complete  solution  cannot  be  obtained  in  closed  form. 

sour* 

dn 

dw 

lo  -w/"i 

“l 

V1  " 

■y2w1/m  . 

We  find 

n(0) 

V1 

(119) 

£ * 

D 

1 Vl 

v 167T3/2e2r0  ’ 

(120) 

E(0) 

■ Vl6Tr3^2mr0v1/3  , 

(121) 

0 * 

y/mr0vl/3YT  , 

(122) 

Ns  - 

J JJ  ■ ^ n(0,t°  ■ 

(123) 

As  a function  of  the  dimensionless  potential 

♦ * 

-e$(i) 

w ’ 

1 

(124) 

the  number  density  and  electric  field  are 

n “ 

n(0)[e*^  - Ytt  erfc  (VTJO  ] 

(125) 

E * 

E(0) [ (1  - 2ip)e-4'  ♦ 2 VtT4>'>/2  erfc  (Vlp)] 1/2  . 

(126) 

The  inverse  of  iKz) 

is 
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♦ 

z “ "^1  *0  ^/" t ( 1 - 2t|/)e-^  + ZVttV^2  erfc CV1JT) ] ~1^2  (127) 
0 

Here,  erfc(x)  is  the  complementary  error  function  1 - erf(x).  Figures  3, 
4,  and  5 show  n(z)/n(0),  E(z)/E(0),  and  \)<(z)  obtained  from  a numerical 
evaluation  of  Equations  125  through  127.  The  abscissae  in  Figures  3,  4, 
and  5 are  scaled  to  the  Debye  lengths  Equations  74,  90,  105,  or  120, 
appropriate  to  each  example. 


The  asymptotic  behaviors  are 

2£2 


; 4 oo  n 

n(z)— n(0)  -f  , 

(128) 

z 

7 - « 

B(*)- -E(O)  , 

(129) 

2 

- — r , 

* ft  * 

(130) 
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so  that  4>  diverges  logarithmically.  For  z <<  we  find 


n(z)  - n(0)  - 1.602^‘J  , 

E(z)  * ECOJ  [l  - VT^]  • 

t- 


(131) 

(132) 

(133) 


so  that  dn(z)/dz  diverges  at  the  surface. 


The  dipole  moment,  Equation  62,  diverges  logarithmically  with  *t> 
and  is  only  weakly  dependent  on  fluence,  as  was  the  case  in  Example  3. 
P(z)  vs.  z is  shown  in  Figure  6. 
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Figure  6.  Dipole  moment  vs.  distance  from  surface. 


An  electron  ejected  with  z component  of  velocity  vzQ  and  wz(0)  ■ j 
‘q  has  z component  of  kinetic  energy  at  z given  by 


wz(z)  = wz(0)  + e$(z)  , 


(134) 


which  is  shown  in  Figure  7 for  several  values  of  w (0),  showing  how  various 

z 

particles  slow  down  with  distance  from  the  surface.  The  maximum  distance  such 

a particle  will  go,  zmnv(w  (0)),  and  the  time  t (w  (0))  to  reach  this 

ni3x  z max  z 

distance  are  shown  in  Figure  8.  Particles  ejected  with  wz  < Wj  go  less  than 
two  Debye  lengths. 


The  total  stored  kinetic  energy  is 
K = 3 K 


= 1.012  w N , 
1 s 


(135) 


where  K„  is  the  total  kinetic  energy  stored  in  the  z component  of  motion, 
and  the  energy  in  the  electric  field,  Ug,  is  2K/3  in  accordance  with  Section  6. 
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Figure  8.  z^x  and  *max  as  a function  of  normal 
emission  energy. 


K(o)) 


4. 


Values  of  the  Debye  lengths,  Equation  120,  are  given  in  Reference 


Referring  to  Figure  3,  we  see  that  in  this  example  the  number 
density  drops  off  very  rapidly  at  first,  falling  to  half  its  surface  value 
in  less  than  1/4  Debye  length.  Both  of  the  finite  energy  spectra.  Examples 
1 and  2,  give  a very  poor  representation  of  the  electron  density.  In 
Example  3 the  absence  of  the  low-energy  electrons  causes  the  density  to 
drop  off  more  slowly  than  in  Example  4. 

The  low-energy  electrons  double  the  surface  value  of  n.  Equations 
104  and  119.  The  electric  fields  have  a very  similar  behavior  as  seen  in 
curves  3 and  4 in  Figure  4. 

In  Example  4,  Figure  6 shows  that  the  dipole  moment  varies  slowly 
with  distance,  remaining  of  order  Wj/4ne  for  many  Debye  lengths.  But  at  one 
Debye  length  it  is  only  0.11  Wj/47ie. 


SECTION  8 


SUMMARY 


Ne  have  presented  the  solution  of  the  steady-state  one-dimensional 
boundary  layer  problem,  summarizing  previous  work  and  presenting  the  solution 
for  the  interesting  exponential  energy  spectrum  in  some  detail. 

Tn  general  we  find  both  the  electron  number  density  at  the  surface 
and  the  surface  charge  density  (and  also,  therefore,  the  electric  field  at  the 
surface)  are  independent  of  the  electron  charge.  The  characteristic  potential 
of  the  layer  and  the  dipole  moment  per  unit  area  depend  only  weakly  on  the 
fluence  in  the  case  of  an  exponential  energy  spectrum.  General  principles 
imply  that  the  energy  stored  in  the  electric  field  is  just  twice  that  stored 
in  the  kinetic  energy  of  the  normal  component  of  motion. 

Estimates  for  the  times  at  which  the  steady-state  solution  should 

apply  show  it  to  be  valid  if  the  X-ray  flux  changes  only  slightly  in  about 

-3  2 

1 nanosec  for  fluences  above  10  cal/cm  . 
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